The existence of a bending rigidity for a hard sphere hquid 
near a curved hard wall: Helfrich or Hadwiger? 
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^ • Abstract 
J? . 

^ In the context of Rosenfeld's Fundamental Measure Theory, we show that the bending 

rigidity is not equal to zero for a hard-sphere fluid in contact with a curved hard wall. 

f-— ^ , The implication is that the Hadwiger Theorem does not hold in this case and the 

^ ' surface free energy is given by the Helfrich expansion instead. The value obtained for 

, the bending rigidity is (1) an order of magnitude smaller than the bending constant 
g ■ 

I I associated with Gaussian curvature, (2) changes sign as a function of the fluid volume 
! 

Ch ■ fraction, (3) is independent of the choice for the location of the hard wall. 

O 
O 



I. INTRODUCTION 

The Helfrich free energy |l| has proven to be an invaluable starting point in the 
description of the surface properties of complex surfaces such as membranes or sur- 
factant systems It is the most general form for the surface (or excess) free 
energy of an isotropic surface expanded to second order in the surface's curvature: 



X 

= jdA[a - 5a J + kK + ..], (1) 

where J =1/ Ri + 1/ R2 is the total curvature, K = l/ {R1R2) is the Gaussian curvature 
and i?2 are the principal radii of curvature at a certain point on the surface. 
The expansion defines four curvature coefficients: a, the surface tension of the planar 
interface, 5, the Tolman length [3], k, the bending rigidity, and k, the rigidity constant 
associated with Gaussian curvature. The original expression proposed by Helfrich 
[l| features the radius of spontaneous curvature Rq as the linear curvature term 



(Sa 2k /Rq 



but in honour of Tolman, who was the first to consider curvature 



corrections to the surface tension 



J], we use the notation in Eq.([T| 



ace the Helfrich free energy in certain 



si based on the imphcations of the 



Recently, an alternative description to rep 
situations was put forward by Konig et al. 

Hadwiger Theorem jo], [lo|. The Hadwiger Theorem states that any functional of a 
system that is translationally invariant, additive and continuous, can be written as 
a linear combination of the four Minkowski functionals: volume, surface area, and 
the integrated total and Gaussian curvatures 10|. The implication is that, as an 



alternative to Eq.([T]), the surface free energy can be written as: 

f^Hadwiger = JdA [ a - 5a J + k K] . (2) 

Comparing the two expressions for the free energy in Eqs.([T]) and (j2]), we are led to 
the following two implications of the Hadwiger Theorem: 

1. The bending rigidity constant is zero, 

2. Higher order curvature terms, represented by the dots in Eq.([T]), are absent. 

The question now is for which systems are the conditions of the Hadwiger Theorem 
fulfilled so that the bending rigidity and higher order curvature terms are all strictly 
zero? It was suggested that for a hard sphere fluid in contact with a hard, structure- 
less wall, the Hadwiger 



its surface free ener 



?heorem should hold and Eq.© is a complete expression for 
irgy 13, M . The evidence for this suggestion is based on a numer- 

n n n 

ical analysis p, |8|, Oj of the free energy in spherical and cylindrical geometry using 
Rosenfeld's Fundamental Measure Theory (FMT) 12, fsj, showing that the bending 
rigidity is zero within numerical accuracy. To understand the basis for this result in 
more detail, we consider surfaces for which the curvatures J and K are constant. The 
Helfrich free energy per unit area is then given by: 

nu/A = a{J,K) = a-5aJ+^J^ + kK + ... (3) 

For a spherically or cylindrically shaped surface with radius R, this expansion then 
takes the form: 

25a {2k + k) / X /.X 

^s{R) = ^ - ^ + ^2 + • • • (sPlie^e) (4) 

a,{R) = ^ - ^ + ^ + • • • (cylinder) (5) 

Note that only the combination 2k+k appears in the expression for the surface tension 
of the spherical interface, so that the conclusion whether k is identically zero or not 



can be made only from an analysis of the cylindrical system. Next to the curvature 
dependent surface tension, one may also investigate the curvature dependence of the 
wall density . According to the wall theorem, the wall density of a fluid in contact 
with an infinitely hard, planar wall is related to the bulk pressure p through an ideal 



gas law [14 



(6) 



kBTp" =p 

The wall theorem can be generalized to a spherically shaped hard wall jlS^lTj or to 
a more generally shaped hard wall with (constant) curvatures J and K [18i] : 

k 



kBTp'^{J,K) = p + aJ-25crK--r + 2kJK + ... 

Note that a term proportional to is absent in the expression above [1 
spherically or cylindrically shaped surface, this expansion takes the form: 



(7) 
For a 



kBTpfiR) 



2a 26a 



k^Tp-^iR) = p + 



a 



k 



+ ... 



(sphere) 
(cylinder) 



(9) 



R 2R^ 

where the dots represent terms of 0{1/R^) which indicates that the term proportional 
to 1 /i?^ is absent in the expansion of the spherical interface. The corresponding term 
in the expansion of the cylindrical interface is related to the bending rigidity thus 
supplying a second route to the determination of its value. Note that these expressions 



are valid only when the radius R is defined via the wal 
In this article, we revisit the analysis by Konig et al. 



density p^ = p{r = R'^). 
71, [Sj for a hard sphere fluid in 

7i,y,y, 



contact with a hard wall. Using the exact same theoretical model as in refs. 
,.e. FMT HQ. we sKow ,n SecfoMthat a detaUed .,.e*a. analysis y.elds a 

bending rigidity that is not equal to zero, but an order of magnitude smaller than 
the rigidity constant associated with Gaussian curvature. Consistent values for k are 
obtained from the analysis of the radius dependence of the surface tension, Eq.([5]), 
as well as from the analysis of the radius dependence of the wall density, Eq.([9]). 

As a further consistency test, we perform a systematic expansion of the FMT free 
energy to second order in the curvature for the spherical and cylindrical interface in 
Section IIIII T his expansion is analogous to a similar expansion for the liquid- vapour 
interface |20[ . It is shown that the resulting expressions for a, S, and k are all in terms 
of the fluid density profile of the planar interface, po{z), whereas the expression for 
the bending rigidity k, features the leading order curvature correction to the density 



profile, pi{z). The values obtained for a, Sa, and the combination 2k + k using these 
expressions are all consistent with the results of Konig et al. j?, 8| and those by Bryk 



et al. 



11| . The value obtained for the bending rigidity is not zero and consistent with 



the two values obtained from the radius dependent surface tension and wall density, 
urthermore, it is in qualitative agreement with recent MD simulations by Laird et al. 



2l\ who determined the curvature dependent surface tension of a fluid near a hard 



wall by Gibbs-Cahn integration 



21 



22|. 



II. FUNDAMENTAL MEASURE THEORY 



In this section, we discuss Rosenfeld's Fundamental Measure Theory [12] as it is 
applied specifically to a one-component fluid consisting of spherical particles with a 
diameter d. The free energy is then the following functional of the fluid density p(r) 



12 



13|: 



knT 



df 



P ln(p) - P 



u Vcy-t (r) 



(10) 



kBT k^T 

where /x is the chemical potential and where the external field V^xt(^) is used to express 
the presence of the hard wall. For spherically shaped fluid particles the free energy 
density = 0(n2, ns, n^) is explicitly given by 

1 



-722 ln(l - rig) H — ^ h 



2(1-^3) 



24(1 -n3)2 



The three densities naif) (a = 2, 3, V) are different convolutions of the fluid density 



na{ri) = df2 p(f2) Wain - 



(12) 



where the weight functions Wa{r) are exphcitly given by |l3[ | 



d d T d 

W2{r) = S{- - r) , W3{r) = e{--r) and wv(r) = - 5(- - r) . (13) 

The Euler-Lagrange equation that minimizes the free energy in Eq.( ITO|) is given by 



^ -ln(p) ■ ^-^^'^ 



(14) 



ksT kBT dna{f2) 

Note that the Euler-Lagrange equation features Wa{r2 — ri) and not Wa{ri — r2) as 
in Eq.([nD (l3|. 

For a uniform system, we have that n2 = 6r]/d, ns = rj and ny = 0, with the 
volume fraction defined as rj = (tc/G) pd^. The Euler-Lagrange equation in Eq.([Ti 



then becomes: 

^ ,n , , r/ (14 - 13r/ + 5r/^) 
= - ^"(^ - ^) + 2(1^^^? • 

Using the expression for the chemical potential above the bulk pressure is obtained 

from Q = —pV giving the Percus-Yevick equation of state: 

k-^T 7r(l — f])'^ 

We mention that a refinement of FMT was recently proposed 23| to yield the more 
accurate Carnahan Starling equation of state ^] instead of Eq.f lT6|) . It is expected 
that results do not depend sensitively on this refinement. 

Next, we consider the implementation of FMT in three different geometries: the 
planar, spherical, and cylindrical interface. 

Planar interface 

In planar geometry, we can simplify the expressions for na{r) = na{z), where z is the 
coordinate normal to the interface, as: 

riaizi) = jdz2 p{z2) wli^zi - Z2) , (17) 
where the weight functions w'l{z) are explicitly given by 

w',{z) = vrrfe(^-l^l), (18) 
wliz) = vr(|-.^)e(^-|.|), 
w^z) = 2nze{^-\z\), 

and where nylr) = nY{z)z. The Euler-Lagrange equation in Eq. fll4l) simplifies in 
planar geometry to 

The external field mimics the presence of a hard wall for z<0, i.e. Vext{z) = 00 when 
z < and zero otherwise, so that the density p{z) = for z<0. The surface tension 
is the surface free energy per unit area {a = {Q + pV)/A [25.] ) : 



00 

cr 



dz 



kBT 

-d/2 



(20) 



5 



where the lower integration reflects the fact that <^ and the convoluted densities na{z) 
are zero only when z<—d/2. 



Spherical interface 

In spherical geometry, the densities na{r) = na{r) , with r the radial distance, are: 

^2(n) = Jdr2 p{r2) w^in - r2) , (21) 
n3(ri) = Jdr2 (^^^ p{r2) w'^{ri - , 

nv(ri) = Jdr2 (^^^ p(r2) Wy{ri - rg) + ^ Jdr2 (^^^ p(r2) w'^^ri - rs) , 
where the weight functions are equal to those in planar geometry (Eq.( fT8|l ): 

w',{n-r2) = ndei^-\n-r2\), (22) 

d^ d 
w^in - ra) = 7r(— - (n - rz)^) 0(- - In - ra]) , 

Wy{ri - rs) = 27r(ri - rs) <d{^ - |ri - rg]) , 
and where ny{f) = nY{r)r. The Euler-Lagrange equation in Eq.( lT4l) now reduces to 

_^.M,),l|^,S/*.(^)5^<(.-.) (23) 

Again, the external field mimics the presence of a hard wall, i.e. V"ext('") = oo when 
r < R, which serves to define the location of the radius R of the spherically shaped 



hard wall. The surface tension now becomes: 



oo 



k^T J \RJ 

R-d/2 



dr — 



: 



k^T k^T 



(24) 



Cylindrical interface 

In cylindrical geometry, the densities ^a(r') = nQ,(r), with r the radial distance to the 
cylinder axis, reduces to: 

?22(ri) = Jdr2(^^y p{r2)w2iri,r2) , (25) 
nsin) = jdr2(^^ p{r2)w'^{ri,r2), 



nv(ri) = Jdr2(y-^ p{r2) Wy{ri,r2) + ^ /^^^ (^) '^(^2) ^s'l^b ^2) 



where the weight functions are given by: 

w',in,r2) = 2dK{f3)Q{^~\r2-n\), (26) 

wl{n,r,) = 16rir2[i?(/3) + (/32-l)K(/3)]e(^-|r2-ri|), 

w'y{n,r2) = 16nr2[K{/3)-E{/3)]e{^-\r2-n\), 

Wyin, ra) = 4(ri - r2)K{l3) e(^ - |r2 - ri|) , 

where n\{r) = nY{r)f and where K{I3) and E{P) are complete eUiptic integrals of 



the first and second kind, respectively [26|. The argument of the elliptic functions 
is defined as = — (r2 — Ti)^]/(4rir2). Note that the weight functions in the 
cylindrical case are functions of the radial distances ri and r2, separately and not 
only the difference ri — r2. The Euler-Lagrange equation in Eg. (1141) in cylindrical 
geometry reduces to 

Again, the external field mimics the presence of a hard wall for r < R. The surface 
tension in cylindrical geometry is given by: 



oo 



ksT J \R 

R-d/2 



p\nip)-p-/-p + <l>+/-Qir-R) 



{2i 



The procedure to evaluate (Ts{R) and (Jc{R) is now as follows. For a certain fixed 
value of the fiuid volume fraction r], the corresponding chemical potential and pressure 
are determined from Eqs. OlSp and (1161) . Next, a value for the radius R is chosen and 
the Euler-Lagrange equation in Eq. fl2^ or fl27|) is solved numerically to obtain the 



density profile p(r) (for details, see the excellent review on FMT by Roth in 13|]). The 
density profile thus obtained then directly provides the wall density p^ = p{r = R^) 
and the radius dependent surface tension by evaluating the integral in Eq. fl2^ or 
(1281) . Finally, the curvature coefficients are obtained from a fit of the surface tension 
and wall density plotted as a function of and comparing with the expansion in 
Eqs-dl]) and (E]) or Eqs.(|H]) and (Q. The fit is carried out by varying the reciprocal 
radius from to 0.1 in steps of 0.01 in units of 1/d. The resulting 11 data points are 
then fitted (least-square) to polynomials in 1/R of progressing order starting from a 





a 


Sa 


2k + k 


V 


Eq.dMD 




Eq.dMD 




Eq.dMl) 


0.10 


-0.0220978 


-0.00130941 


-0.00130941 


0.000428974 


0.000428971 


0.20 


-0.1394516 


-0.014812 


-0.014811 


-0.0014205 


-0.0014207 


0.30 


-0.512482 


-0.07321 


-0.07318 


-0.0161 


-0.0161 



TABLE I: Numerical values for the surface tension a (in units of k^T/cP), Tolman length 
6(7 (in units of k^T/d) and the combination 2k + k (in units of k^T) for three values of the 
volume fraction 7]. These vahies are determined from an analysis of the radius dependence 
of the surface tension of a spherical interface and by a direct evaluation of the expression 
in Eas. (p3]l - (l35]l . 

quadratic polynomial to a polynomial of order 7. It is verified that the coefficients in 
the fit level off with the variation used as an indication of the numerical error. 

For the spherical interface, the polynomial fit of crs{R) provides values for the 
coefficients a, 6a, and the combination 2k + k. The results are listed for three fluid 
volume fractions in Table [11 The values for a and 6(7 obtained from the polynomial 
fit of the wall density are, within error, equal to those listed in the Table. 

For the cylindrical interface, the polynomial fit of (Jc{R) again provides values 
for the coefficients a and 5a (which are consistent with the results in Table [T]), but 
the coefficient of the l/i?^-term now yields values for the rigidity constant k. These 
values are not equal to zero within numerical accuracy and are listed separately 
in Table [Tll Also listed are the values obtained from the polynomial fit of the wall 
density. Already it is noted that these two approaches are consistent and lead to the 
conclusion that the bending rigidity is not equal to zero for this system. To further 
corroborate this result, we consider a third approach in the next section. 

III. CURVATURE EXPANSION 

In this section, we expand the free energy of the spherical and cylindrical surface 
systematically to second order in 1/R. The analysis is outlined explicitly for the 
spherical interface - the analysis of the cylindrical interface is more or less analogous, 
but we indicate where it differs from that of the sphere. 



8 





bending rigidity k determined via: 


V 






Eq.dMD 


0.10 


0.000415172 


0.000415165 


0.000415171 


0.20 


0.00074260 


0.00074251 


0.00074254 


0.30 


-0.000685 


-0.000615 


-0.000619 



TABLE II: Numerical values for the bending rigidity k (in units of k-^T) for three values 
of the volume fraction rj. The bending rigidity is determined in three different ways: by an 
analysis of the radius dependence of the surface tension of a cylindrical interface, the radius 
dependence of the fluid wall density of a cylindrical interface, and by a direct evaluation of 
the expression in Eq. (f39l) . 

Spherical interface 

All quantities are expanded to second order in the curvature. In particular, the 
expansion of the density Ps(r) reads: 

Ps{r) = po{z) + + + ... , (29) 

where z = r — R. The coefficients in the curvature expansion of the density are 
determined from the curvature expansion of the Euler-Lagrange equation in Eq.( !23ll . 
The result is that the (planar) density profile po{z) is determined from Eq. flTU]) and 
Pi{z) follows from solving: 

= ^ + E h^ f,of^]t\f , n\{zr)wl{z,-z,) (30) 
Po(^i) t^-^ a<(2;2)an0(z2) ^ 

where 0o = 0({'^a}) and where we have defined z^ = zi — Zx. As we show below, it 
turns out that for the evaluation of the curvature coefficients it is sufficient to obtain 
the density profiles Pq{z) and p\{z^ only. Using the expanded density profile, we can 
then determine the coefficients in the expansion of n^(r): 

<(r) = <(.) + ^ + ^ + ..., (31) 
where ^^(2;) is given by Eq. (fT7I) and n\[^z) can be calculated from 

^2(2^1) = p2;2 [pi(2;2) + ^12Po(^2)] ^2(^1 - ^2) , (32) 



'^v(^l) = jdZ2 [[Pl{z2) + Z12 Pq{z2)] wl{zi - Z2) + Po{z2) wl{Zi - Z2)'\ . 

Again, the evaluation of n^(2;) tuns out not to be necessary. 

The expansions for Ps(r) and ^^(r) are inserted into the expression for the surface 
tension in Eq. (!24l) . Making a systematic expansion to second order in 1/-R, using the 
Euler-Lagrange equations in Eqs.( fT9|) and ( l30l) . one uhimately obtains expressions 
for the curvature coefficients by comparing to the curvature expansion in Eq.(jll). For 
the surface tension of the planar interface the result in Eq. (!20|) is recovered: 
a 



dz 

d/2 



Li p 

Po In(po) - Po - -j-^ Po + 00 + -j-Tf Q{z) 



For the Tolman length one obtains 
6a 



knT 



dz z 

d/2 



Po In(po) - Po 



^ Po + 00 + T^efzl 



knT 



knT 



(33) 



(34) 



00 00 



-d/2 

For the combination 2k + k one finds 

— rvn 

2k + k 



" -d/2 ^ 

1 °° °° d(f> 

"2 l^'^hdAW)''^"'^""'^"''"'^- 



kuT 



dz z^ 

d/2 



p p 

Po In(po) - Po - Po + 00 + -j-Tf 
fceT kBT 



(35) 



oo oo 



+ V fdzi fdz2 Zl ° . Po{z2) zuwKzi - Z2) 

" -d/2 ^ 

4e J-.,....,..:..-.. 



-d/2 
00 00 



n]^{zi) po{z2) Z12 w^{zi - Z2) 



^ 00 00 
+ 2 jdzi jdz2 



-d/2 
00 00 



■pl{z2)wl{Zi - Z2) 



+ jdzi jdz2 

-d/2 



dn%{zi) 



Po{z2) ^2^3(^1 - Z2) . 



10 



By solving the density profile po{z) from Eq. flT^ and pi{z) from Eq. fl5U]) . these 
coefficients can all be evaluated directly without having to determine the full radius 
dependent surface tension as a function of 1/R. It is therefore no surprise that this 
route to the evaluation of the curvature coefficients is much more convenient. To 
compare our results to the results by Bryk et al. listed in their Table I jll|, we need 
to take care of the fact that in their analysis the location of the radius is defined 
according to the location of the "actual surface" which accounts for the fact that the 
molecule's center of mass is half a diameter away from the surface when it interacts 
with the hard wall, -Ractuai = R — d/2. The curvature coefficients are then shifted 
according to the following transformation: 

(^)R-d/2 = ^ + ^ ' (36) 
,5. , 5- pd'^ ad 



'K-d/2 - - g 2 ' 

(2k + k) = 2k + k + — + 6a d . 



The form of these transformations are derived by shifting the location of the z = 
plane by a distance d/2 in the expressions in Eqs. fl33l) - fl35l) . 

The results for a, 6a, and the combination 2k + k are plotted in Figure [T] as the 
solid lines. Also shown in Figure [1] are the calculations from Bryk et al. [11] (circular 
symbols), computer simulation results by Laird et al. 2l|, |22| (square symbols) and 
Scaled Particle Theory (SPT) j^^] (dashed lines), for which the expressions read: 

ad^ 3ri{2 + ri) 



k^T 27r(l-?7)2 
6a d Si] 



knT An (1 - r/) ' 

2* + * -±,„(i-„) 



(SPT) (37) 



ksT 47r 

From the results in Figured] it is concluded that the curvature coefficients calculated 



using Eqs. fl33|) - fl35|) are consistent with those obtained by Bryk et al. ll|], although 
there seems to be some small discrepancy for the combination 2k + k at larger volume 
fractions. We come back to this point in the Discussion. 

Cylindrical interface 

The analysis for the cylindrical interface is more or less analogous to that of the spher- 
ical interface, with one notable difference being that the weight functions w'^{ri,r2) 

11 



0.5 





0.2 ^ 0.3 




0.2 ^ 0.3 




0.2 ^ 0.3 



FIG. 1: Various curvature coefficients as a function of tlie fiuid volume fraction: (a) surface 
tension a (in units of k^T/cP), (b) Tolman length 6a (in units of k^T/d) and (c) the 
combination 2k + k (in units of k^T). The drawn lines are the results calculated using 
Eqs. (p3|) - (j35|) . transformed according to Eg . ([36]) so that the radius is defined as that of the 



"actual surface" . Circular symbols are results from Bryk et al 
computer simulation results by Laird et al. 



21 



11[; square symbols are the 



22l |: the dashed line is the SPT result in 



12 



in Eq.( 12B]) also need to be expanded in 1/R: 



W2(ri,r2) 


= W2izi - 


^2) 1 


/ 




= ^3(^1 - 


^2) 1 




W3,(ri,r2) 


= ^3(^1 - 


^2) 1 








■^2) 





d^A - z 
16i?2 

3(^74 - 4) 



'^ + ...'1, (3^ 



32i?2 

2 

12 



+ 



^2/4-^2 



+ 



16 R'^ 

Following the same procedure as for the spherical interface, the expressions for cr 
and 6a in Eqs.f l33|) and are recovered, and one obtains as an expression for the 
bending rigidity k: 

k 



1 f f 

-4 II y y^^2 g-^f^ Pl(^2) Zi2 W^izi - Z2) (39) 
+ i E 7^-1 /rf-2 g n(,^t°(..) "^("^^ - 



dn^{zi)dn'li{zi] 



-d/2 

+ 1 E Jdz, /d.2 g n.,^t°(..) - 



-d/2 
00 00 



+ i E /rfz, /d.2 P0(.2) - H2) - Z2) 

" -d/2 ^ 

-1 Jdz, /d.2 UZ2) (^ - ZL) -3°(-l - -2) 

-d/2 ^' 

-d/2 ^ 

1 °° °° (90 

+2 hhd^) 



p0[Z2) 2:12^3(^1 - ^2) 



-d/2 

where pi(2;) and nl^{z) are the same as in the spherical analysis. It is noteworthy 
that since no reference to the location of the z = plane is made in this expression, 
the bending rigidity is independent of the choice for the location of the radius R, i.e. 
(A;)R_d/2 = k. In this respect the bending rigidity is a much more inherent property 
of the interface in question. The result of the evaluation of the bending rigidity 
using Eq. (1391) is shown as the solid line in Figure O The open circles and crosses 
are the previous results for k listed in Table [Tll Also shown are very recent computer 
simulation results by Laird et al. [21!] (solid circles). 
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0.002- 



0.001- 



k 0- 



-0.001- 



-0.002- 




0.2 ^ 0.3 



FIG. 2: Bending rigidity k (in units of k^T) as a function of the fluid volume fraction. The 
drawn line is the result calculated from Eq. (j39p : open circles and crosses are the previous 
results listed in Table [III solid circles are the computer simulation results by Laird et al. 
21|. 



Figure [2] is the main result of this article. It shows that the bending rigidity is 
definitively not equal to zero in the context of FMT theory for a hard sphere 
fluid near a hard wall and the Hadwiger Theorem does not apply in this case. We 
have shown this via three more or less independent approaches which agree within 
numerical accuracy with each other. A further corroboration of this result are the 
computer simulation results by Laird et al. 2l|]; although the agreement is not 
quantitative, the shape of the volume fraction dependence of k is strikingly similar. 

Finally, we like to mention that by combining the expressions in Eqs.( |35ll and ( |39ll . 
an expression for the rigidity constant associated with Gaussian curvature may be 
obtained: 

k 



knT 



dz 

d/2 



po In(po) - Po - Po + 00 + 



knT 



kuT 



(40) 



oo oo 

Jdzi jdz2 Zi 



-d/2 



dnl{zi) 



Po(^2) Zi2W^^{Zi - Z2) 
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oo oo 



+ Jdzi jdz2 zi ^ po{z2) wl{zi - Z2) 



-dj2 

00 00 



- J2 Jdzi jdz2 (°^) /'o(^2) - 34) wl{zi - Z2) 



-d/2 
00 00 



1 /" /" dcj^Q dj^ 

-d/2 

Note that k can be evaluated from the properties of the planar interface only; a result 



that is consistent with similar expressions for the liquid-vapour interface 



20|. 



IV. DISCUSSION 

We have shown that the bending rigidity is not equal to zero in the context of FMT 
theory for a hard sphere fluid near a hard wall and that the Hadwiger Theorem does 
not apply in this case. Evidence for this conclusion is shown in Figure [2] where the 
results of three independent approaches are shown to agree within numerical accuracy. 
Noteworthy is that the bending rigidity changes sign from positive to negative as a 
function of increasing fluid volume fraction. It is smaller than the rigidity constant 
associated with Gaussian rigidity roughly by an order of magnitude. 

The reduced magnitude of k may certainly be partly responsible for the fact that 
in a previous analysis j3] it was hard to distinguish it from zero. Another possible 
source for the discrepancy may be due to a different fit procedure used to extract 
the curvature coefficients from the radius dependence of the surface tension and wall 



density. A comparison between our analysis and the analysis in refs. [?, 8, 11 1 shows 
that while numerical results for as{R) agree to within a high degree of accuracy Q], 
the difference in fit procedure leads to a fitted value for 2k + k that may differ by as 
much as 10 % (see Figure [T]:). One may very well speculate that the difference in fit 
procedure used may also have consequences for the fitted value obtained for k. 

The question now remains, what is the underlying physics of the Hadwiger Theo- 
rem? The Hadwiger Theorem is not merely some abstract notion from Mathematics 
and one should be able to understand more microscopically when the conditions (i.e 
additivity) that lead to it are fulfilled. To address this question, let us consider 
the general form of the mean-field expressions for the surface tension in spherical 
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(Eq.(^) and cylindrical geometry (Eq.i^) |29j: 



'^siR) = ldr{^Ju,{r), (41) 



a,{R) = Jdr n,(r) 



where n(r) depends on the distribution of the fluid density p(r) in the interfacial 



region and may be referrec 
the excess lateral pressure 



to as the excess free energy density or (the negative of) 



29| 



Now, if it is assumed that the lateral pressure is independent oi R, i.e. n(r) = 110(2;), 
then the only radius dependence in Eq. fHT]) is due to the geometric factors (r/i?)^ 
and (r/R). Therefore, we immediately conclude from Eq.( HT|) that 



a = Jdz Uoiz) , (42) 
6a = — Jdz zIlo{z) , 
k = Jdz z"^ no(^) , 

and the bending rigidity k is zero. Furthermore, all the higher order terms in the 
expansion in 1/i? are absent. [It was already Helfrich himself who derived these "ge- 
ometrical expressions" in terms of progressing moments of the excess lateral pressure 
[3^.] These results correspond precisely to the predictions of the Hadwiger Theorem 
so that we may conclude that the Hadwiger Theorem corresponds to the statement: 

Hadwiger Theorem <^=^ H(r) = Ho (-2) . (43) 

This means that the Hadwiger Theorem applies when the fluid molecules do not 
rearrange themselves when the curvature of the interface is changed. Certainly, the- 
oretical models may be constructed in which such a rearrangement does not occur, 
but in general this is certainly not the case. To explore this curvature dependence, 
we expand Hs(r) for a spherical interface in 1/R: 

U,ir) = Uo{z) + ^ + .... (44) 

Szleifer and coworkers jsi] already showed that the bending rigidity k is then ex- 
pressed as 

k = ^JdzzUi{z), (45) 

which explicitly demonstrates the conclusion that k results from the (possible) rear- 
rangement of molecules when the curvature of the interface is changed. An example 
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FIG. 3: Density profiles po{z) (in units of and pi{z) (in units of as a function of 

z (in units of d) for 77 = 0.8. The values at z = correspond to the pressure po{0~^)=p/kBT 
and (twice) the surface tension pi{0~^) = 2a / k^T; cf. Eq.®. 

of such a rearrangement of molecules as described by the density profile pi (z) is shown 
in Figure [3] for 77 = 0.3. 

Now, one could argue that the vanishing of the bending rigidity is simply a matter 
of length-scale j?]. The length-scale associated with the molecular rearrangement due 
to curvature is the width of the interfacial region ^, which is small compared any 

n 

to macroscopic length-scale unless the system is critical p^] or when a macroscopic 



wetting layer is present 



32|-I34j| . However, the same argument would apply to all 
the curvature coefficients and in particular to the rigidity constant associated with 
Gaussian curvature which scales similarly to the bending rigidity. 
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